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Abstract 

Casimir effect, in a broad interpretation which we adopt here, consists in 
a backreaction of a quantum system to adiabatically changing external con- 
ditions. Ahhough the system is usuaUy taken to be a quantum field, we show 
that this restriction rather blurs than helps to clarify the statement of the 
problem. We discuss the problem from the point of view of algebraic structure 
of quantum theory, which is most appropriate in this context. The system 
in question may be any quantum system, among others both finite as infi- 
nite dimensional canonical systems are allowed. A simple finite-dimensional 
model is discussed. We identify precisely the source of difficulties and infini- 
ties in most of traditional treatments of the problem for infinite dimensional 
systems (such as quantum fields), which is incompatibility of algebras of ob- 
servables or their representations. We formulate conditions on model ideal- 
izations which are acceptable for the discussion of the adiabatic backreaction 
problem. In the case of quantum field models in that class we find that the 
normal ordered energy density is a well defined distribution, yielding global 
energy in the limit of a unit test function. Although we see the "zero point" 
expressions as inappropriate, we show how they can arise in the quantum 
field theory context as a result of uncontrollable manipulations. 
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1 Introduction 



The Casimir effect, bearing its name from the pioneer work by Casimir has 
become in recent decades an increasingly popular topic in quantum field theory, 
with a new review of the subject appearing every few years, see |21 El EI El El • The 
effect consists in the response of a quantum field, even in a ground state, to the 
introduction of external, usually macroscopic, bodies. Initially the effect existed 
as a theoretical prediction only, and a rather mysterious one, for that matter. 
However, increasing experimental evidence of its existence (see e.g. has lead 
to attempts at better understanding of its theoretical foundation. The problem is, 
that the theoretical side of the phenomenon has been plagued from the beginning 
by divergent expressions, as well as conceptual difficulties, which have proved to be 
surprisingly persistent. This is the more surprising, that models usually considered 
in this context are linear, so the usual sources of quantum field infinities are absent 
here. 

In an earlier paper |7I I have given a diagnosis of the reasons of this state of 
affairs and proposed to treat the problem from the algebraic point of view. This is 
the most natural and fruitful framework in quantum physics, with its beginnings 
already in the classical book on quantum mechanics by Dirac, and modern devel- 
opments in quantum field theory and statistical physics described e.g. in mono- 
graphs |H| and in] . When viewed from that angle the source of difficulties is rather 
obvious, and can be briefly termed as uncritical use of the concept of quantum 
field 0. More precisely, what we mean is this. The first step to define a quantum 
theory is to identify a set of quantum observables (we ignore here the question 
of non-observable variables) together with algebraic relations between them, such 
as canonical commutation relations. Once we have this, a concrete physical real- 
ization of the theory corresponds to a choice of a representation of the algebra 
of observables. Non-comparable physical situations are realized by non-equivalent 
representations p . Although we want to see the real world as a unity, physics, 
of course, is about idealizations, and various idealizations need not be compatible 
(take e.g. an isolated system and a thermodynamic limit system). However, if we 
want to consider transitions from one physical situation to another and compare 
values of one and the same observable in various states, all situations taken into 
account must be describable in one common representation. Now, these scheme 
is violated in most treatments of the Casimir effect. For a typical situation of 
a quantum field in a region with movable sharp boundaries the difficulty arises 
already on the algebraic level: there is no consistent choice of an algebra of ob- 
servables for all physical situations coming into play. The energy of the "free" 
field is an observable defined in the vacuum representation of the algebra of the 
field smeared with Schwartz test functions. For this algebra evolutions imposed 
on the field by the presence of boundaries cannot be defined. Furthermore, even if 
one "smooths out" the boundaries so as to make a common choice of an algebra 
possible, one still has to satisfy rather severe restrictions necessary to ensure the 
equivalence of representations. These restrictions are typically violated in usual 
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treatments. For these reasons we have advocated in [7] the view, that the model of 
sharp boundaries, as well as many other insufficiently regular models, are wrong 
idealizations in the context of Casimir effect, and we have also proposed (and an- 
alyzed) a class of models imitating Dirichlet conditions. Let us stress this: once 
a model has correctly been chosen, there is no space (nor need) for further ad hoc 
regularizations, and the formalism should yield well-defined answers to legitimate 
questions. Although views on nonphysical nature of sharp boundary conditions 
have been also expressed elsewhere (see e.g. ^01)) it seems that the conditions for 
a model to be acceptable in the sense described above have not been analyzed 
before. For instance, in a series of recent papers Graham et al. 'llj investigate 
a linear model imitating Dirichlet conditions. Being linear, the model should be 
well-defined without any renormalization (except for a trivial normal ordering for 
quadratic quantities like energy density; for external potentials without bounded 
states this is our example (iii)i at the end of Section 13 below). However, renor- 
malization is ad hoc imposed on it by the authors in order to give meaning to 
a meaningless expression. 

The algebraic problems we have described often do not appear if one restricts 
attention to local quantities in quantum field theory. This fact is connected with 
what in algebraic formulation is called the local quasiequivalence of representations 
(see |S|). The point is as follows. In quantum field theory observables are equipped 
with the property of locality: each local observable carries as a label an open space- 
time region with compact closure in which it may be measured. As stated above, 
two representations of the totality of these observables representing two different 
physical situations may be non-equivalent. However, physically one would expect 
that even if the two situations are globally non-comparable, one should be able to 
compare results of local measurements (think of the vacuum representation and 
a thermodynamic limit representation). Mathematical formulation of this expec- 
tation is this: think of states in each of the representations as density operators; 
restrict attention to an arbitrarily chosen compact region of spacetime; then for 
each state in one of the representations there is a state in the other which yields 
the same expectation values for observables localized in the chosen region. If the 
two representations have this property they are called locally quasiequivalent. As it 
turns out physically important representations do indeed often have this relative 
property, and then expectation values of local quantities may be compared. How- 
ever, in the situation we want to consider in this paper this result falls short of 
our needs in twofold way. First, we want to calculate expectation values of global 
quantities, which are limits of local ones for the size of the spacetime region tend- 
ing to infinity - in this case the global differences of the representations come into 
play. Second, in situations like fields with imposed boundary conditions, even finite 
regions which overlap with boundaries are not local in the above sense: for those 
regions even the scopes of local algebras in presence of boundaries are different 
than in the vacuum theory. 

Another important point we want to stress in our analysis of Casimir effect 
is the choice of the observable to be compared in various considered states. In our 
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view the backreaction of a system perturbed by external agents is determined by 
the expectation value of the energy as defined by the unperturbed system, one 
and the same (as an operator) in all states to be considered. A more systematic 
discussion of this point in a wider context will be found in the next section. Here 
we want to note that some local, in the spirit of the last paragraph, calculations 
of the Casimir energy do follow similar ideology; in the gravitational context see 
esp. a paper by Kay jjl], and for electromagnetic field with conducting boundaries 
a paper by Scharf and Wreszinski |13| . However, in many other local calculations, 
esp. those using "the Green function method" , the situation is somewhat ambigu- 
ous: it is often not clear enough what the general viewpoint is, and the result may 
agree with the above method in some cases, but disagree in others. We shall dis- 
cuss this point more fully in Section below. For the global energy, as determined 
by the unperturbed system, to be defined in states of the system influenced by ex- 
ternal conditions, as required by the above ideology, we need one common algebra 
and globally equivalent representations, as explained earlier. This imposes restric- 
tions on the perturbed dynamics, which are usually violated, and the transition 
from local Casimir energy to global one is then blocked by infinities of physical 
nature. Any "regularization" thereof is an ad hoc procedure, striving at this late 
stage to compensate for the wrong idealization in interaction with external con- 
ditions. Finally, there is a group of works explicitly comparing the expectation 
values of different global observables: energy with and without interaction. Here, 
apparently, is the place of the "zero point" ideology. We shall come back to this 
point later on, here we only note, that in this case infinities are even more likely 
to appear. In that method one subtracts expectation values of regularized "bare" 
energy observables; different energy observables may have different singularities, 
not cancelling under subtraction. 

The present paper is the first of the two in which we develop and describe 
more fully what was announced in [7] (we use notation slightly changed at some 
points with respect to that paper). Here we discuss more general results on the 
admissibility of models for the purpose of investigation of quantum backreaction. 
In the second paper applications to particular models are discussed. We use rig- 
orous mathematics, and present real proofs. However, we hope that the paper is 
readable for a wide audience. 

In Section 121 we place the quantum field Casimir effects in a wider context of 
a backreaction of a quantum system to adiabatic changes. This section thoroughly 
discusses the foundation for the calculation of this backreaction in any quantum 
system. In Section |21 we discuss quantization of a class of linear systems, which 
include quantum fields under linear external perturbations. We put stress on less 
widely known aspects of this otherwise standard procedure which are important 
in the present context. Section ^ discusses an application to a finite-dimensional 
system. In Section |S1 we treat infinite-dimensional cases, and we formulate condi- 
tions for admissibility of a model for the discussion of backreaction effects. More 
specifically, we consider a quantum field case in Section O We show that with 
a slight strengthening of these conditions not only global energy, but also energy 
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density may be defined, and in the appropriate limit global energy is recovered. 
Section Ul contains somewhat more explicit discussion of the points made earlier in 
this Introduction on the existing calculations of Casimir effect. We also comment 
there on the "zero point" expressions for the Casimir energy. We try to under- 
stand, from the point of view of the formalism presented in the present work, how 
such expressions may arise. We show how imposing unacceptable idealization of 
sharp boundaries and doing unjustified manipulations leads from our expression 
for the energy density to "zero point" expressions for Casimir energy. Appendix 
gives a simple form to a handful of mathematical facts in Fock space which are 
needed in the main text. These are known results, but we believe that this sum- 
mary makes some of them more accessible. 

2 A quantum system under external conditions 

Trying to put the discussed phenomenon in a broader context we shall adopt the 
following point of view. The Casimir-type effect consists in the backreaction of 
a quantum system on the adiabatically changing external conditions under which 
the system is placed. 

The background for this idea is this. We consider a larger closed system con- 
sisting of two subsystems Q and M. These subsystems interact with each other, 
but to certain degree (this will be made more precise below) maintain their sepa- 
rate identity. Part Q is our relatively simple quantum system under consideration 
(say, electromagnetic field) , while M is supposed to be of much more complicated 
nature (say, conductor plates), and to have among its variables some of collective, 
macroscopic type (separation of the plates). We want to determine the effect of 
the evolution of the joint system on the collective variables attached to M . 

Because of the complicated nature of the part M of the system and its in- 
teraction with Q, to tackle the problem one has to make some simplifying as- 
sumptions. There are at least two possibilities, both of them of phenomenological 
nature. In both cases one simply represents part M of the system by a few col- 
lective variables (such as separation of the plates), suppressing all the details of 
this subsystem, and representing the interaction between M and Q by some sim- 
ple effective model. The first possibility is to equip the collective variables with 
a fully quantum nature, and put forward a simple model for the closed system. 
This approach, when applied to the more specific situation of a quantum field 
in interaction with macroscopic bodies, is chronologically more recent one in this 
field, and is called the dynamic Casimir effect (see 0). Although we admit that 
this forms an open possibility, we shall not take it up in this article. Firstly, not 
much can be said with high degree of certainty and mathematical rigor. Secondly, 
the apparent attractiveness of the approach does not necessarily withstand a closer 
scrutiny. A macroscopic body undergoes "constant observation" , so effects of de- 
coherence play primary role, which is not taken into account in this approach. 
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Another possibility, which we take up in this paper, has more restricted as- 
pirations, but admits mathematically rigorous results, as we are going to argue 
below. We have to admit, however, that there is some confusion at its physical 
formulation. We hope to contribute to its removal. This second approach consists 
in approximating the collective quantities, which characterize a macroscopic body 
as a whole, by classical variables. Moreover, one considers only situations, in which 
the whole system changes adiabatically. The effect of the evolution on the macro- 
scopic (classical) variables in this context is what we referred to as a Casimir-type 
effect at the beginning of this section. More specifically, the Casimir effect refers 
to a quantum field in interaction with macroscopic bodies. 

One should be more specific about physical assumptions and approximations 
involved in the situation implied in the last paragraph. This is, in our opinion, 
a point not clear enough in many discussions of the Casimir effect. Therefore we 
shall try to be systematic, even at a risk of being too detailed. 

(i) One considers first the isolated quantum system Q {M is absent). We give 
its description in the algebraic formulation of the Heisenberg picture, see 
e.g. 0. 

(i)i Basic quantum variables at a fixed time form an abstract *-algebra A, 
e.g. an algebra of canonical commutation relations (or, more technically, 
its exponentiation to the Weyl form). 

(1)2 This algebra is represented by operators in a Hilbert space H: 

7::A^n{A), A^7:{A), n{A*) = 7r{A)\ 
TT{aA + pB) = aTi{A) + I3ti{B) , tt{AB) = ii{A)t:{B) , ^ ' ^ 

where 7r(^) is a concrete algebra of operators in Ti.. Vectors in that 
space, or, more generally, density operators acting in this space, rep- 
resent states of the system Q. Representation tt is assumed to be irre- 
ducible; then vectors correspond to pure states. 
(1)3 The intrinsic dynamics of Q is defined by an automorphism of the al- 
gebra A: 

at-.A^A, A^atA. (2.2) 

This automorphism is implemented by a unitary evolution in the Hilbert 
space H: 

n{atA) = U{t)T:{A)U{ty , U{t) = e^p{itH) , (2.3) 

where H has the interpretation of the energy operator of the system. 
This operator is supposed to have nonnegative spectrum, and usually is 
assumed to have a ground state, represented by a unit eigenvector to the 
lowest point in the spectrum. One does not perturb the above relations 
by adding a multiple of the identity operator to H, so the ground state 
may be assumed to have zero energy. By irreducibility of tt the energy 
operator H is then uniquely determined. 
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(ii) One introduces now part M into the system. This part is characterized by 
classical variables (we shall denote them by a), so no new quantum vari- 
ables are added. Therefore system Q should retain its identity, and changes 
in its state will influence the classical variables of M . Thus various states 
to be considered must be physically comparable. These assumptions have 
mathematical consequences. 

(ii)i Identity of the system Q is formed by the algebra A ((i)i above), so 
this algebra must remain unaffected by M. 

(11)2 Physical comparability of states demands that also the particular rep- 
resentation TT of ^ ((1)2 above) remains unaffected by the introduction 
of M. 

We stress the importance of this point as it is both crucial for the scheme, 
as wc sec it, and usually overlooked. If various physical situations to be 
considered demanded different algebras or different (nonequivalent) repre- 
sentations, the approximation would break down, as one could not follow 
the change in the system Q brought about by the creation of (and changes 
in) M, and its reaction to that occurrence. Further support for this point will 
be found below. Let us note again, what was discussed in introduction, that 
the local quasiequivalence of representations if Q is a quantum field system 
is not enough for our purposes. 

(iii) We consider now dynamics in presence of M, and assume at first that the 
variables a are frozen. In this case Q is still a closed system in interaction 
with conditions created by M, and for each fixed a its evolution is again 
given by an automorphism of the algebra A: 

aat--A^ A, aatA . (2.4) 

One assumes implement ability of new evolutions in the representation n: for 
each a we have 

n{aatA) = Ua{t)7r{A)Ua{t)* , C/„(t) = eM^Ha) • (2.5) 

For each a the generator Ha is defined by this up to the addition of a multiple 
of the identity operator, so we have the freedom 

Ha^Ha + Xaid, (2.6) 

where Xa is any real function of parameters a. 

(iv) One allows now the coupled system Q — M to evolve. Part Q alone is not 
a closed system any more, so it could be too restrictive to assume that the 
evolution of its variables would be given at the algebraic level, as an automor- 
phism. However, this evolution should still be describable in terms of unitary 
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operators in the Hilbert space Ti. (not forming a one-parameter group, in gen- 
eral); this corresponds to the assumption of conservation of probabiHties in 
the subsystem Q. The use of the Schrodinger picture for the quantum part Q 
wiU be more convenient in the present context. State of the coupled system 
Q — M is specified at a given time by a vector in the Hilbert space H. (de- 
scribing the state of Q), and values of a and, possibly, their time derivatives. 
We formulate the evolution of this system. 

(iv)i Suppose that a{t) is known as a function of time. We assume that this 
functional dependence is very slow (system M is "heavy"). It is then 
a justified approximation to assume that the time-dependent hamilto- 
nian of the evolution of the system Q is given by i?a(t) (with Ha defined 
in (iii) above). As a{t) is slowly varying we assume the adiabatic ap- 
proximation to calculate the evolution. One is usually interested in the 
situations in which the initial state of Q is given by an eigenvector of Ha 
for the initial value of a. Suppose that for each a we have a nondegen- 
erate, normalized eigenvector ipa of Ha'- 

Ha^a = Ea^a , (2.7) 

and the family ipa depends continuously on a. If at f = the state of Q 
was given by V'a(o): then at later times in the adiabatic approximation 
its state is equal to ip{t) — e*'^'*Va(t): where (p{t) is a real function 
depending functionally on Ea and ipa- If an operator B represents an 
observable, then the time-dependence of its expectation value is given 

by 

{B)t^{Mt),B^ait)), (2.8) 

so it is a function of a in this approximation. It is important to note 
that the eigenvalues Ea are modified by the addition of Xa under the 
transformation l|2.6|l , but both the eigenvectors ipa and the mean values 
{B)t remain unchanged. 

(iv)2 Finally, the evolution of the macroscopic variables a{t) must be de- 
termined. This is the most controversial part of the problem, but we 
believe that the foregoing discussion indicates its proper solution. 
The intrinsic energy stored in the quantum part Q is represented (in the 
Schrodinger picture) by the operator H {{1)3 above), which in the cou- 
pled system is not a constant of motion any more. Under the assump- 
tions of (iv)i its expectation value is a function of a, depending on the 
choice of the continuous family of eigenvectors ipa- 

£a:= i^a,H^a), (2.9) 

and the time-dependence of this expectation value is through a{t) only. 
Changes in £a correspond to the energy which has been transferred 
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from Q to the rest of the system, which (with the suppression of all 
microscopic details of M) is described by the variables a. Thus £a plays 
the role of a potential energy with respect to these variables. We assume 
that the rest of the total energy of the coupled system is supplied by 
the kinetic energy of M, thus we obtain a potential system, with the 
generalized force given by 

df 

Ta = —^. (2.10) 

With a specific form of the kinetic energy for a particular model the 
motion of a(t) could be determined, and with large inertial parameters 
(a "heavy" system) the approximation of its slow change should be 
confirmed. 

We have thus spelled out all the assumptions and arrived at the basic formu- 
las (|2.9() , (|2.10() . In the following sections we shall take these formulas as a starting 
point. The derivation of the formulas was not rigorous, as this would demand more 
information on the underlying microscopic model of the closed system Q—M. A de- 
tailed analysis of these questions is both outside the usual discussions of Casimir 
effect, and also beyond the reach of a rigorous calculation at present. However, 
we believe that the proposed discussion offers more plausibility than most of the 
statements of the problem to be found in literature. In particular, points made 
by us in (ii) above are typically ignored; we shall see their consequences when Q 
is an infinite-dimensional system, e.g. a quantum field. Furthermore, we want to 
draw a closer attention to the formula l|2.9|) and contrast it with what one obtains 
by the generalization of the "zero point" method to the more general context dis- 
cussed in this section. In the latter case our formulas (|2.9() and (|2.10() are replaced 
respectively by 

f-P- = - i?0 , K-^ =-^. (2.11) 

where Ea is the eigenvalue determined by H2.7|l , and Eq some reference eigenvalue 
of H . One can object to these formulas on several grounds. 

(a) The philosophy behind them seems to be this: the backreaction of Q on M is 
due to the changes in Ha, which may be interpreted as the sum of intrinsic 
energy H oi Q and some interaction energy. However, we think that it is 
M which absorbs the interaction and transforms it in a phenomenological 
way into an effect on macroscopic variables a, while Q has a rather clear-cut 
identity. 

(b) The energy given by the "zero point" philosophy is not a quantum mechan- 
ical average of any clear-cut observable: with changing a one changes the 
observable Ha- Moreover, as already pointed out. Ha and their eigenval- 
ues are subject to the gauge freedom (I2.t)|l . The usual argument runs that 
this is fixed by the quantization of the "proper" classical expression for Ha- 
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We regard this argument as very unreliable. Quantum theory is the more 
fundamental one, so in case of doubt it should not seek a verdict from the 
classical theory. 

(c) We put forward the following "consistency check" . Suppose that for certain 
values of parameters a the effect of M on Q vanishes. In this case the backre- 
action force should vanish as well. The supposition means that for a = oq the 
vector il)ao is also an eigenvector of H , Hipag = E^ao with some eigenvalue 
E. Using this equation one easily shows that our formulas yield 

■rao—~Cj — U, 

so they pass the check. On the other hand 



a—ao 

which, in general, has no reason to vanish. 

(d) In Section0]below we discuss an example of a Casimir-type effect in a canon- 
ical system with finite degrees of freedom. In this example the "zero point" 
method fails dramatically, yielding a completely unphysical result. 

How, then, may "zero point" expressions arise? We shall show in Section [7| below 
how for quantum fields problems "zero point" expressions may be related to ours 
by unjustified idealizations and manipulations. 
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3 A class of quasi-free systems 

We discuss in this section a general quantization scheme for a class of simple 
models. This class includes linear perturbations of multi-dimensional harmonic 
oscillators or quantum fields. 

Consider first the classical case. Let 7?. be a real Hilbert space, and denote its 
scalar product by (.,.). Let hhe a. selfadjoint, strictly positive (hence invertible, 
with densely defined inverse h"^) operator in 7?,, with the domain Dfi^h). We form 
the external direct sum C = 'Dn{h) (B TZ C TZ (B Tl, and denote its elements by 
V = V (B u, V £ I?7^(/i), u £TZ. With the symplectic form a defined by 

a(yi,V2)^{v2,ui)~ivi,U2) (3.1) 

space C becomes the phase space of a classical model. Let the Hamiltonian function 
of the model be given by Ti{v, u) = ^[(u, u) -f {hv, hv)] (where all mass parameters 
have been absorbed by momenta) . The evolution determined by this Hamiltonian 
in C is given by 

Tt{v Q)u) = { cos{ht)v + sm{ht)h^'^u) ( - sm{ht)hv + cos{ht)u) . (3.2) 
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The differential form of this evolution is actually valid only on a subspace of Z 
(dense in 7?. ® 7?,), but the evolution itself is properly defined on the whole of L. 
Operators Tt form a one-parameter group of symplectic transformations 

TtT, = rt+, , a(TtT/i, TtFs) = o{Vx,V^) . (3.3) 

Note, also, that 

T_f = (id ® - id) Tt (id ® - id) . (3.4) 
Each V' £ C may be identified with an element of the dual space by the rule 

V'{V) = {v\ u) + {u', v) = a{V', (id ® - id)^) . (3.5) 
Then using Eqs. H3.3|l and (|3.4|l one easily shows that 

{TtV'){V)^V\TtV). (3.6) 

The above model may be generalized by considering a more general subspace 
contained in 'Dn{h) © TZ and invariant under the evolution law H3.2|l . We use this 
freedom to choose 

C^Vn{h)®Vn{h-^'^) (3.7) 

(the invariance under H3.2|l is easily checked). The evolution law Tt may be now 
expressed as a unitary evolution in a complex Hilbert space. One introduces a com- 
plex Hilbert space K, which is the complexification oiTZ^ K, — TZ® iTZ, with scalar 
product (denoted by the same symbol) and complex conjugation defined by 

{vi + iUi,V2 + iU2) = (wi, V2) + {ui,U2) + i{vi,U2) - i{ui,V2) , (3.8) 
X ^ Kx = X £ IC , K{v + iu) — v — iu . (3.9) 

We shall write v — Re(f + m), u ~ lm{v -\- iu). The operator h has a unique ex- 
tension to a complex-linear operator on /C, denoted by the same symbol, with the 
domain T>{h) = T>n{h) iDn{h). This new h is again a selfadjoint, positive oper- 
ator, and it commutes with the conjugation. Consider now a real-linear operator 

j : £ Ranj C /C, j{V) = h^/'^v ~ ih-^^'^u . (3.10) 

Its range Ran j is a real-linear subspace of /C, dense in /C, and j is a bijection of C 
onto Ranj. Then for all V € C: 

j{TtV)^e"^'jiV), (3.11) 

so Ran J is invariant under e*'**, and the evolution may be expressed as 

TtV = r'ie'''*jiV)). (3.12) 

Space JC, regarded as a real vector space, has a natural symplectic struc- 
ture introduced with the symplectic form Im(/, 5). Space Ranj' is its symplectic 
subspace. One easily shows that 

a{VuV2)^lm{j{Vi),j{V2)), (3.13) 

so j is a symplectic transformation of C onto Ranj. 
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The mapping j, as well known, serves to construct the ground state rep- 
resentation of the quantum version of the model, and the space K. is then the 
"one-particle space" (see below). A natural problem thus arises: to extend the 
construction of the space C to the largest possible space compatible with the sym- 
plectic mapping (|3.1U|I . that is to extend £ and j so as for Ran j to cover the whole 
space K. (instead of being only dense in /C, as above). One defines on 'Dn(h^^/'^) 
the scalar products 

{vi,V2)+ ^ {h^'^v^,h^'^v2) , vi,V2eVn{h^'^), (3.14) 

{ui,U2)- = {h-^'^ui,h-^'^U2), ui,U2eVn{h-^/^), (3.15) 

and denotes by TZ+ and TZ^ the Hilbert spaces obtained by the completion of 
'D-jiih^^^) and I?7j(/i^-^/^), respectively, with respect to the norms ||u||+ = ^ {v, v)+ 
and = ■\/ (it, u)_ . For V G Vn{h^/^) and u G X>k(/i"^/^) we have 

\\h^/^v\\ = ||u||+ and \\h^^^u\\ = Therefore operators /i^/^ and h'^^"^ ex- 

tend by continuity to bijective isometric operators /i^/^ and h^^/^ respectively. 



h^/^-.n+^n, \\h^/'^v\\ = \\v\\+, (3.16) 

h~^:TZ-^TZ, \\lP^^u\\ = \\u\\^. (3.17) 

We note for future use that 

n±nn^Vn{h^^^^) ■ (3.i8) 



This is easily seen in the spectral representation of h: if ft, is a multiplication by 
a positive, different from zero almost everywhere, function / in a space L^{M, dp), 
then TZ consists of functions i/' for which /^,^ \'\li{m)\^ dii{m) < oo, TZ± consists of 
functions for which Jj^j{f{m))'^^\il){m)\'^dfi{m) < oo, and I?7^(ft,^^/^) - of those 
satisfying both conditions. 

For V € Vn{h^''^) and u G Vnih"^/"^) one has 

\{v,u)\ = \{h^'\,h-^'^u)\<\\vU\\u\\^, (3.19) 

thus {v, u) extends to a continuous pairing 

7^+ X 9 G M, < + . (3.20) 

Now one can set 

£ = 7^+®7^_, ^(1/1,^2) = («2,ui) - (i;i,U2), (3.21) 

3:C^IC, 3(V)^h^^v-ih-^u, (3.22) 

ftV ^r\e^'''3iV)) . (3.23) 
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As a consequence of (|3.18|) one has 



£ n (7^ e 7^) = VTz{h^''^) e Vn{h-^''^) ■ (3.24) 

It is easy to see that now Ran j — IC, the space given by Eq. I|3.7|l is dense in £ (in 
its Hilbert space structure norm), and the time evolution on £ is the continuous 
extension of the evolution on the space 1)3.7(1 . Moreover, 

GiVi)3{V2)) - (1^1,1^2)+ + {u,,U2)-+td{Vi,V2), (3.25) 

so, in particular, j is a symplectic mapping of (£,(?) onto (/C,Im(., .)). Relations 
H3.5|l and (|3.6|l are also generalized to 

V'{V)^{v',u) + {u',v) =ct(F', [id®(-id)]F), (3.26) 
{ftV'){V) =V'iftV). (3.27) 

Once we have the largest arena consistent with the scheme, particular models are 
defined by choosing a subspace invariant under the evolution: 

CcC, ftCcC. (3.28) 

The maximal model is invariant under the time reversal, represented by the oper- 
ator id®(— id) appearing in ((3.26() . We want to retain this property for the model 
defined by C, which is equivalent to the assumption 

C = C+®C-, C±(in±. (3.29) 

Examples of particular spaces include the class of spaces 

C = Vn{h-+^) n Vn{h-^) © Vnihn n Vn{h~'---) , 
r, s, i e (0, 00) , s<t+i, t<s + |, 

all of which are contained in H3.7() . 

The quantum version of the maximal model (with the symplectic space C) 
is now obtained by standard procedure (see e.g. |2, vol. II). Starting with expres- 
sion H3.2(j|l one aims at replacing V' by some "quantum variable" In quantum 
theory a concrete representation of a quantum variable is an operator in a Hilbert 
space. If the classical variable is real, its quantum counterpart should be rep- 
resented by a selfadjoint operator. Thus one assumes that a Hilbert space H is 
given, and for each V G C one has a selfadjoint operator ^{V) in that space. The 
functional dependence of ^{V) on V is assumed to be linear, and the canonical 
commutation relations are imposed: 

mVi),MV2)] ^ ia{Vi,V2)id, (3.31) 
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where one still has to clarify the domain problems, li V = v (B u, then we shall 
also write ^{V) = ^{v,u). The element P{v) = $(w,0) has the interpretation 
of the quantum momentum for the "test vector" v, and X{u) = $(0, u) - of the 
quantum position variable for the "test vector" u. With the linearity of ^{V) the 
above commutation relations are equivalent to those in a more familiar form 

[X{ui),X{u2)] = , [P{vi),P{v2)] = , [P{v),X{u)] = -i{v, u) id . (3.32) 

It is well-known that there are many different concrete representations of the 
above scheme, and this is why it is desirable to formulate the canonical commu- 
tation relations in an algebraic way. As there are no bounded operators satisfying 
these relations, it is usual to take them in an exponentiated variant. This leads to 
the Weyl form of these relations. The Weyl algebra over the symplectic space C is 
the unique C*-algebra generated by elements Vl^(y), V & C, and a unit element 1, 
by the relations 

WiVi)WiV2) - 6-5^(^1' ^2) wiVi + V2) , (3 
WiVy ^W{-V) , T4^(0) = l. 

One looks for representations of this algebra by bounded operators in a Hilbert 
space (which exist for all C* -algebras). Let tt be such a representation in the 
Hilbert space H. One says that this representation is regular, if for each V E C 
the one-parameter group of unitary operators 

R3 t:{W{sV)) (3.34) 

is strongly continuous. If this is the case, then there exist, by Stone's theorem 
(e.g. 531), selfadjoint operators ^{V) such that 

n(W{V)) = exp(i$(F)) . (3.35) 

Moreover, one shows that for each finite-dimensional subspace C C C there exists 
a dense subspace T>' C Ti which is contained in the domains of all operators ^{V), 
V £ C , is an invariant subspace and an essential domain of selfadjointness for all 
of them, and on which linearity of ^{V) in its argument V G C and commutation 
relations H3.31|l are satisfied (this follows from the Stone-von Neumann uniqueness 
theorem, cf. vol.11). While not all canonical systems with these properties arise 
in this way from regular representations of the corresponding Weyl algebra, most 
of those needed in physics do, and one usually restricts attention to this class. 

The algebra A of the maximal model is thus the Weyl algebra over £. Dynam- 
ics of the model is a "quasi-free" evolution obtained by a simple "quantization" of 
the classical evolution Tf. Being guided by the replacement V' ^ ^ in Eq. H3.27|l 
and the relation ()3.35f) . one defines it on the algebraic level by 

at{W{V)) = WiftV). (3.36) 
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One looks now for a representations tt of the algebra in which this evolution law 
may be implemented: 

7:{W{ftV)) = Uit)7riW{V))U{tr , U{t) = e'^^ , (3.37) 

where is a selfadjoint operator. The ground state representation is obtained 
if iJ is a nonnegative operator with zero energy ground state. This represen- 
tation is constructed in standard way with the use of the Fock space method. 
Let Wo(/) = exp[i$o(/)], / € /C, be the Weyl system of operators in the Fock 
space TC built on the "one-particle" space JC (see Appendix IA.1|I . We set 

7r{W{V)) ■.= WoOiV)), (3.38) 

or, which is equivalent, 

7r{W{V)) = e^*(^) , $(y) = M3iV)) . (3.39) 

Using identity 1)3. 25|) and properties of the operators VFo(/) one easily shows that 
this indeed constitutes a representation of the Weyl algebra (|3.33l) . By the ir- 
reducibility of the Weyl system in Fock space this representation is irreducible. 
Moreover, using Eqs. (|3.38(l and H3.23|l one rewrites the condition H3.37|l as 

M^o(e**V) = e**^l^o(/)e-**^, feH. (3.40) 

The discussion of Appendix I A. II shows now that 

H = dT{h), (3.41) 

where dT{h) is the "second quantization" of h (see Eqs. (IA.19l ~ lA.21ll ). This energy 
operator has nonnegative spectrum, and a unique ground state represented by the 
"Fock vacuum" S7. 

Consider now a restriction of this model defined by a subspace C invariant 
under evolution (Eq. H3.28|l ) and time reflection (Eq. H3.29|) ). The algebra of the 
model is the subalgebra of the Weyl algebra (|3.33|) obtained by restricting the test 
vectors to C. It is well-known, that the resulting model is not identical with the 
maximal one \i C ^ C (see jH], vol.11), but one can demand that its ground state 
representation approximates that of the maximal model. This representation may 
be constructed as before, but the scope of Weyl operators used in this represen- 
tation is restricted to {Wo(/) | / £ j('C)} ~ cf. Eq. (|3.38|) . This set is irreducible 
in /C if, and only if, the space j{C) is dense in /C, or, what is the same, C is dense 
in TZ+ ® TZ-. With the time reflection symmetry assumption (|3.29() this takes the 
form 

is dense in TZ+ , £_ is dense in TZ- . (3.42) 

We restrict attention to those spaces C which satisfy this condition. This restriction 
can be paraphrased by saying that there are no superselection rules in the Fock 
space of the ground state representation. 
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Examples 



(i) Multidimensional harmonic oscillator 

In this case 7?. is a finite-dimensional Euclidean space, and his a. positive self- 
adjoint operator defined on the whole of TZ. We choose C — TZ®TZ. Space JC 
is the unitary space obtained by complexification of 7?., and Ranj' = /C. The 
more familiar simple form of the model is obtained by choosing in TZ an or- 
thonormal basis (ei, . . . , e„) of eigenvectors of h and putting Xi — $(0, e^), 
Pi — $(ei,0). The system is then the set of n independent harmonic os- 
cillators with canonical variables {Xi,Pj}, unit masses, and frequencies uji, 
where hci = tOiei. 

(ii) Free scalar field 

Free quantum fields are usually defined as operator-valued distributions on 
test functions of all spacetime variables. The evolution equation (Klein- 
Gordon for the scalar field) is already encoded in this formulation, which 
is manifestly relativistically covariant. For our purposes the equivalent ini- 
tial value formulation is preferable - we want to separate evolution law, as far 
as it is possible, from setting up of the algebra. 

Standard identifications for this model are as follows: 



where subscript M denotes the real part of the respective function space, 
I?(M'^) is the space of infinitely differentiable complex functions of compact 
support, and A is the Laplace operator. Standard solution of the initial value 
problem for the Klein-Gordon equation has the form of Eq. (|3.2|) , and the 
assumptions (I3.28|) . H3.29|l and (|3.42|) are satisfied. 

(iii) Scalar field with external time-independent interaction 

Loosely speaking, the choice of h here is the square root of a selfadjoint 
positive operator of heuristic form "ft,^ = — A + interaction". There are a few 
possibilities. 

(iii)i If the interaction is given by an external field a — a{x) then the choice of 
spaces TZ and K, remains the same as in the free case, while = — A + ct 
(we assume that is still positive - there are no bound states). 
Depending on the form of a the choice of £ as in the free case may 
be admissible (satisfy assumptions (|3.28|) and (|3.42|) ') or not. A safe 
choice for £ is supplied by any of the cases given by Eq. H3.30|l . More 
generally, may be any positive selfadjoint perturbation of —A in the 
sense of operators or forms. 



7^ = L|(]R^), IC^L^iR^), 
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(iii)2 The next possibility arises from restricting the region accessible to the 
field to a proper subset A c K^. In this case TZ = L\{K), K, = i^(A), 
and = —As, where is a selfadjoint extension of the Laplace 
operator defined on twice differentiable functions with support inside A, 
determined by some boundary conditions "i?" . Here, of course, the free 
field choice of C is not admissible, and a safe choice is again given by 
the formula 

(iii)3 Finally, we consider a setting usually assumed for the Casimir effect. 
The whole physical space M'^ is divided by two-dimensional surfaces into 
disjoint open regions Ai, . . . , A^. Position of the dividing boundaries is 
characterized by a set of parameters a. One chooses TZ and K, as in the 
free field case. Depending on parameters a, a family of positive operators 
ha is given by — — A^, where A^, is the Laplace operator in L^(IR'^) 
determined by the assumed boundary conditions (Dirichlet, Neumann, 
etc.) at the dividing surfaces with positions given by parameters a. 
In consequence, the choice of the symplectic space £ must be adjusted 
to the position of the boundaries. A choice of simple possibilities is again 
given by Eq. H^130|l (with ha replacing h). In the Casimir problem one 
wants to compare states of the system at different values of a. However, 
spaces C depend nontrivially on a, thus the respective Weyl algebras 
are also different, and do not define the same quantum system. This 
constitutes the difficulty of traditional treatments of the Casimir effect 
which we anticipated in Section|21 Infinities naturally appear then, and 
are a consequence of an uncritical use of the notion of a quantum field. 

In the following sections we discuss Casimir effects for some systems in the category 
described in the present section. We start with a simple finite-dimensional case. 

4 Deformation of a finite-dimensional harmonic 
oscillator 

Our unperturbed quantum system Q is here an rt-dimensional quantum oscillator 
described in example (i) of the last section. Thus the algebra of the model is the 
Weyl algebra based on a finite-dimensional symplectic space C = TKBTZ, and its 
representation is given by Tr{W{V)) = Wo{j{V)) in the Fock space Ti based on 
the finite-dimensional one-excitation space /C. The energy operator of the model 
is given hy H = dT{h). 

We consider now a combined system Q — AI, as described in Section 12 and 
assume that the influence of M on Q for frozen parameters a manifests itself in 
the change of axes and frequencies of oscillations. Thus the time evolution of the 
algebra for frozen parameters is given by 

aatiW{V))^W{Tat), (4.1) 
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where Tat has the form H3.2|) , but with operator h replaced by an operator ha from 
a family {ha}. For each a the irreducible representation of the algebra in which this 
evolution is implemented by a unitary one-parameter group, with a nonnegative 
energy operator, is constructed by the same method, as in the free Q case, in the 
same Fock space H. Thus 

Ja{V)^hl/^V-th-'/^U, 7TaiW{V)) = Wo{ja{V)). (4.2) 

The Hamilton operator is given by dT(ha), with ground state described by the 
Fock vacuum il. However, as discussed in Section|21 we want to describe the same 
physical situation with the use of the representation tt. Therefore for each a we 
look for a unitary operator Ua which by similarity transforms representation tt^ 
onto tt: 

Ua7Ta{W{V))U:=TT{WiV)), UaMVP: = "^(V) , VGC. (4.3) 

We substitute here the definitions of the representations tt and tTq, and set 
ja{V) = f. This condition then takes the form 

UaWo{f)U:^Wo{Laf), /e/C, where i„ (4.4) 

Both j and ja are bijective symplectic mappings, so La is a symplectic transfor- 
mation of the space (/C,Im(. , .)), and the above condition states that Ua imple- 
ments the corresponding Bogoliubov transformation in the Fock space (see Ap- 
pendix lX!3|l . As K. is finite-dimensional, such Ua exists. The explicit form of trans- 
formation La is easily obtained: 

Laf^ h^'^ha^'^ Re / -f ih-^/^hl'^ Im / , (4.5) 

and then from ljA.29|) one finds La — Ta + Sa, Ta complex-linear and Sa complex- 
antilinear, 

Ta = \{B-^+B:), Sa^\{B-^~B:)K. (4.6) 

where 

Ba = h]J^h-^/\ Kf^f. (4.7) 
In the representation vr the Hamiltonians of the new evolutions are given by 

Ha = UadT{ha)U: , (4.8) 

and the ground state of Ha is given by 

= Ua^ . (4.9) 

Suppose now, that a is a single real parameter, and under the influence of 
the external conditions the state of the subsystem Q changes adiabatically over 
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the states J7a, as discussed in Section 12 The potential for the backreaction force 
is therefore, in accordance with point (iv) in Section |21 determined by 

£a^ina,Hna)- (4.10) 

We take into account that H — dT(h) and use the expression for a form matrix 
element of dT{h) as given in Eqs. HA.23p . (|A.24p : 

f,=^||a(/ii/2/0f^a|P, 

i 

where {fi} is an arbitrary orthonormal basis of /C. We use Eas. l|A.40p . ljA.42|l 
and (|A.45p to find 

Thus we obtain 

£a = Tr [h^/^SaSa*h^^^] = ^ Tr [{ha - h)h-\ha " h)] . (4.11) 
Let the eigenvalues and orthonormal eigenvectors of h and ha be given by 

hCi — C^Ci , haCai — ^ai^ai ■ (4-12) 

Then using the spectral representation h~^ = e~^\eak){eak \ and employing the 
basis Bi for the calculation of the trace we find 

£,=^l!!^J^|(e.,,e.)p. (4.13) 

i,k 

We consider a simple example. Let Q be a two-dimensional oscillator in phys- 
ical space, vectors 61,62 being its main axes, and let the effect of the external 
conditions be the rotation of these axes by an angle ip (= a) , without a change in 
the frequencies. In this case we have e^pk = e^, (61^1, 62) = —(6,^2, 61) = sinip, so 

£^= ^''~^''^\ e-'+e-')sm'^. (4.14) 



The backreaction "force" in this case is a torque 

\2 



^ (62 — ei) I -\ , . „ 

T^^ = (ei +62 )sm2^. 



(4.15) 



For ei — !■ (with 62 kept constant) the torque tends to infinity. This is what one 
should expect. This limiting case describes the situation in which the harmonic 
force in the direction of 62 extends translationally invariant in the direction of 61; 
any rotation of this picture involves an "infinite" change. 

Note, that the "zero point" prescription for the force gives zero in the above 
example, which is an utterly unphysical prediction. 
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5 An infinite-dimensional system 



Let now TZ be an infinite dimensional real Hilbert space. We want to consider 
a situation analogous to that discussed in the last section: system Q defined by 
the operator h, and its perturbations by a family of operators ha- We need to 
take into account complications arising from the unboundedness of the operators, 
as explained in Section |21 

We want to be able to define for our model both evolutions (that determined 
by h and by ha), and both ground state representations. Thus the model has 
to fit into structures defined in Section |31 both by h as well as ha- In particu- 
lar, its symplectic space should in a canonical way be a part of both C and Ca- 
However, the construction of these spaces is based on different, in general, parts 

of 7^ ® 7^ {Vn{h^''^) ® Vn{h-^''^) and Vn{hl^^) ® Vniha^^^) respectively), and 
without some restrictions there is no canonical way of identification of their parts. 
We assume that 

V± = Vnih'^^^^)r]VTz{ht^^^) is dense in 7^± and in 7^a± , (5.1) 

(in this, and similar statements below, signs are either all upper, or all lower). 
In this case space TZ± is the completion of 'D± with respect to the norm ||.||±, and 
TZa± is the completion of the same subspace with respect to ||.||a±- 

Suppose now that the space of a model contains at least a subspace such 

that 

e £" , £^ C I?± , is dense in 7^± and in 'Jla± , (5.2) 

which is a strengthening of the condition (|3.42|) . Note that then j{C'^) — j{C^), 
ja{CP) — ia{^^), and both spaces are dense in /C. Under these assumptions we 
show that the following three conditions are equivalent: 

(i) The symplectic mapping 

ia:=JJa':ja(/:")^j(/:°) (5.3) 
extends to a bounded operator in /C, with a bounded inverse. 

(ii) The operators h and ha satisfy the conditions 

2?7^(/i±i/2) ^ p^(;,±i/2) ^ (5.4) 

Ba = h]J'^h~^/'^ and extend to bounded operators in JC . (5.5) 

(iii) There exists a selfadjoint, positive, bounded operator Ca in TZ, with bounded 
inverse, and such that 

ha = h^^^Cah^^^ (5.6) 

in the sense of forms, that is: ha is the unique selfadjoint operator defined 
by the closed form q{vi,V2) — {h^^^vi,Cah^^^V2) with the form domain 
Qiq)=Vn{h'/^). 
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If these conditions are satisfied, then Ca = C, so both evolutions are well defined 
in C. 

Note that (i) is a necessary condition for the ground state representations 
defined by h and ha to be equivalent. The equivalence implies that Eq. H4.3|) is 
satisfied in particular for all V € £°, or, equivalently, Eq. (|4.4|l for all / € ja{C°)- 
But then, as shown in the Appendix IA.3I La extends to a bounded symplectic 
mapping in /C, with a bounded inverse. 

Let (i) be satisfied. Then on £° from Eq. (15. 3|) we have Laja — j and 
ja — L^^j, which implies that for w± G Cj_ one has 

||/i±i/2u;±|| < const. \\h^^/^w±\\ , \\h^^/^w±\\ < const. \\h^^/^w±\\ . (5.7) 

This means that the norms ||.||± and ||.||a± are equivalent on C^., so they yield 
the same completion, on which these inequalities are preserved. But is dense 
in TZ± and TZa±, so TZ± — TZa± as sets, hence also £ = £o as sets, with equivalent 
norms. Moreover, from Eq. H3.18|l we find 

'Dn{ht^''') =na±nn = n±nn = Vn{h^^''') ■ 

The boundedness of Ba and Ba^ follows now from (|5.7() . which ends the proof 
of (ii). Conversely, if (ii) is satisfied, then one easily shows that the formula 
Laf = B~^ Re/ + iBl Im/ gives the extension needed in (i) (cf. Eq. (|4.5I) '). 

The equivalence of (ii) and (iii) follows by polar decomposition of closed oper- 
ators (e.g. mi). If we assume (ii), then h]/"^ = Bah^^^ with the domain I?7j(ft,^/^) 
is a selfadjoint operator, so ha ~ h^^^B^Bah^^^, and Ca — B^Ba fulfills the condi- 
tions of (iii). Conversely, let Ca be bounded, positive selfadjoint, with a bounded 
inverse. Then Cl^^h^^^ with the domain equal to V-jilh^/'^) is a closed opera- 
tor. Indeed, let w„ e X>7^(/ii/2), ||t;„ - and ||Cy^/(>/2(i;„ - v„,)\\ -> 0. 
But Ca^^'^ is bounded, so also ||/i^/^(wn — Wm)|| — * 0. As h^^^ is closed, there exists 
V € Vqilh}/'^) such that — v\\ — > 0, ||ft^/^(w„ — u)|| 0, and by boundedness 
of also \\Cl''^h^/'^{vr,-v)\\ 0, which shows that Ca^'^h^^'^ is indeed closed. 
Thus the form q defined in (iii) is closed, and the condition H5.6|l means that 
|Cy^ft.^/^| = h]J'^ . It follows that there exists an orthogonal operator Fa such that 
cU'^h^l'^ — Fahi^^, so hi^^h^^/'^ extends to a bounded operator Ba = F*cV^ ■ 
As Ca has a bounded inverse, so the same is true for Ba^ , which ends the proof 
of equivalence of (i) - (iii) . 

These preliminary results show that if the two ground state representations 
are to be equivalent in our model, we have to assume that (|5.6() . and then all 
conditions (i) - (iii), are true. Then the space C is invariant under both evolutions 
and forms the widest possible space of the model. Any subspace (if it exists) 
C = £+ ffi £_ C C which is also invariant under both evolutions and dense in C 
can also be taken as the symplectic space of the model. 
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With these assumptions the symplectic mapping La decomposes into the 
bounded complex-hnear and complex- antilinear parts, La = Ta + Sa- The apph- 
cation of the results described in Appendix IA.3I shows that the necessary and 
sufficient condition for the unitary equivalence of the ground state representations 
is that Sa is a Hilbert-Schmidt (HS) operator, that is 

Ma = Tr [SaS*a] < oo . (5.8) 

Going through the steps (|4.1UI - I1.11|) in the present infinite-dimensional context 
we see that £a is finite if, and only if, S^h^^'^ extends to a HS operator, and then 

£a = Tr [h'^'SaS*ah'/^] . (5.9) 

The quantity TVa appearing in (|5.8I) has a clear-cut physical meaning. The 
results of the Appendix I A . 31 show that if H5.8|l is satisfied, then Ua G 'D{N), where 
N is the "particle" (excitation) number operator. A calculation analogous to that 
carried out for the energy yields 

{na,N^la)=Na, (5.10) 

SO Ma is the mean value of the excitation number in the ground state. 

In the rest of this section we obtain the following criterion for admissible 
perturbations. Let ha be given by (|5.t)|) . The ground state representations are 
unitarily equivalent {Na < oo) if, and only if, 

a = id -1-5,, (5.11) 

where 5a is any operator satisfying conditions 

5a is a HS operator, id-f^a > . (5-12) 



In this case we can write in the sense of forms 



ha = h + h^'^5ah^''^ ■ (5.13) 
Moreover, if conditions (|5.12|) are satisfied, then £a is finite if, and only if, 

dah^^^ extends to a HS operator. (5-14) 
With the condition H5.12|) satisfied one has 



and if in addition H5.14|) holds, then 



■Tr 



id +5a 



\A+5a 



(5.15) 



(5.16) 
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To prove these assertions note first that equations H4.6|l remain in force with 
our assumptions in the present infinite-dimensional context, and then 

SaS: = \^^^^. (5.17) 

If the ground state representations are equivalent, then S* is a HS operator, so 

Ca^^'^iCa — id) is HS as well. But is bounded, therefore also Sa = Ca — id is 
HS. The second condition in (|5.12l) is satisfied by the positivity of Ca- Conversely, 
suppose that Sa satisfies conditions (|5.12|1 . By the first of these conditions Sa has 
a purely discrete spectrum with no other convergence points than zero, and then 
by the second Ca = id +Sa > b id, with 6 > 0. Hence is bounded, and Ca ^^^Sa 
is HS, so Ma = Tr [SaS:] = i TT[6l{id+6a)'^] < oo. 

If the conditions (|5.12|) are satisfied, then in a completely analogous way one 
proves the equivalence of the condition H5.14|l with the finiteness of £a, and the 
equation (|5.16(l . 



6 Energy density of quantum field 

In this section we consider the case of a quantum field, and for definiteness 
we take the scalar field (massive or massless) with standard commutation rela- 
tions and free evolution. Thus here TZ = i^(]R^), h — \/m?- id —A, (m > 0), 
and we take for the test function space the largest space £, as described in 
the previous section. Perturbations ha are assumed to be in the class defined by 
Eqs. IjS.lll - IO^ (in fact, a slight strengthening of these conditions will be needed). 
We show that in this setting the energy density in the ground states fia is well de- 
fined as a tempered distribution, and for the test function tending to unit function 
one recovers the energy expectation value £a- 

We assume a slight strengthening of our assumptions and demand that for 
some a € (0, 1) there is: 

/^(i+a)/4^^/^(i+a)/4 is a HS operator . (6.1) 

Note that this statement with a = is a consequence of our earlier assumptions. 
Indeed, if Sa and Sah}^"^ are HS, then 

Q<TT{h^'Hah'"f= lim Tr[P<o,„>(M(/»'/''5„/ii/4)2p^o,„^(/j)] 

= Jiim^Tr [Sah^'^Sah^'''P(o,n){h)\ = "^^{Sah^'^f < oo , 

where {Ppih)}, F a Borel set in M, is the spectral family of h. 

Loosely speaking, the energy density operator of the scalar field is determined 
by the point-splitting procedure and normal ordering with respect to the vacuum 
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as 

H{x) - lim : H^ix, y) : = lim [H2{x, y) - (f^, H^iS, y)n)\ , (6-2) 

y^x y^x 

H2{x, y) = ^{P{x)P{y) + VX(f) • VX{y) + m^Xix)X{y)) , (6.3) 

where X{u) = $(0,u), P{v) — $(u,0) (see H3.32|l and the preceding remarks), and 
to get X{x) and P{x) one sets formahy v and u equal to Dhac deha concentrated 
at X. We are interested in the energy density (fla, Hix)fla) in the ground state fla- 
We now make this precise. The real Schwartz test function space Sr is con- 
tained in 2?-fj(/i^/^) n T>Ti(h~^/'^), so functions from that space may be used for 
"smearing" both X{x) as P{x). Let 'Wi,W2 & 5h. The precise meaning of H6.3|) is 



^2(^1, k;2) 

= i($(wi,0)4>(w;2,0) -^^(O, Vwi) • $(0,Vw2) + m2$(0,i(;i)$(0,u;2)) . (6.4) 

To find normal-ordered expectation value (fJa, ■ H2{wi,W2) '■ ^a) we need to know 
(r^o, ■^{Vi)^{V2) : ria), where we assume that G iSr ® Sr. We recall the defini- 
tions of the representations tt and tt^ and their equivalence relations: 

(Eqs. (|4.2(l . 14.9(1 and 14.:-?|l respectively). Using them one finds 

(f},$(i^i)$(F2)r!)-^(j(^i),j(^^2)) 

(6 5) 

= l{h'/'v^,h'^'v2) + \{h~^'^U,,h-^/^U2) + ^^(^1,^2) , 
{na,^{Vi)^{V2)na) = [n,<^>a{V^)^a{V2)n) = \{]a{Vl),]a{V2)) 

(6 6) 

= \{hl'\,,h\'^V2) + l{h-'/'m,h-'/'u2) + '-^{V,,V2) , 

{n,,Mv,)HV2):n,) - ]^{h^/\,^5,h'/\2)~\{h-^'^u,,r^h~''^u2) , 

(6.7) 

so 

{VLa,:H2{wi,W2) ■■ i^a) = T'a(wi, W2) = Tal{wi,W2) + Ta2{wi,W2) , (6.8) 

where 

1 ^2 

Tal{w^.W2) - - f/l^/'W, TX7V/^'/'«^2) , (6.9) 

4 V id +da ' 
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rr I \ ^fl,i/2 ,1/2 \ rn /,_i/2 ,-1/2 \ 

Ta2(wi,W2) ^ -:{h ' wi, „ h ' W2] —{h ' wi, h ' W2 

4 V id +d„. / 4 V id +o„ / 



4 V id +Sa 



(6.10) 



We have added the conjugation sign over wi on the r.h. side to make the expression 
hnear rather than antilincar also for complex functions. The Tai part and the first 
term in Ta2 result from splitting 

+-r^- (6-11) 



id +6a id +Sa 

We show that: 

(i) Ta{wi, W2) defines a distribution Ta{x, y) on S 

(ii) For each 7/ S the expression Tq (^ + ^, ^ — ff) is a distribution on S 
and for each test function / the function 

7? ^ 8a{v, f)^ j Ta{i + rj, f - rf)fiO d'^ (6.12) 

is continuous and bounded (we use the "integral" notation of distributions). 
The energy density according to point-splitting procedure is then the distri- 
bution 

£a{f) = £a{OJ) . (6.13) 

(iii) Let MO = /(eC), /(O) = 1, / G S{R). Then 

\im Saife) ^ £a ■ (6.14) 

Before starting the proof we fix conventions for the Fourier transforms. 

For a,b e M" we set 

fib) - (2^)-"/2 / /(a)e-'^-"d3„ , f{a) = /(-a) . (6.15) 



We consider the Tai and Ta2 parts separately. Expressions H6.12|) and H6.13|l for 
Tai replacing Ta will be denoted £ai{v^ f) ^^'^ ^ai{f) respectively. 

As ^6a{id+Sa)~^^'^h^^^ is a HS operator in L^(]R^), it is an integral operator 
with a kernel ka{x, y) £ i^(R^) (see e.g. [21 )• Thus Tai is obviously a distribution 
on 5(R^), determined by the ordinary function 



Tai{x,y) = I ka{z,x)ka{z,y)d z. (6.16) 
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As for each ff there is ka{z,^ + ff)ka{z,£^ — fj) E {M.^ , cP z , the distribution 
^aiiVi /) is indeed well defined, 

SaMf)^ JkJxI+^)fiOka{z,^-ri)d^zd^^. (6.17) 

Now, Fourier-transforming ka{z, £, + ff) and f{£^)ka(z, £, ~ fj) with respect to z and 
^ one finds 

I ka{r,P)Ur,q)f{p-'f)e-''^' ^P + '^^d^rd'pd^q, (6.18) 

as the integrand on the r.h. side is absolutely integrable. Therefore £ai(?7, /) is 
continuous in ff. For 77 = we get 

£al{f)= l\ka{z,OffiOd'zd'C (6.19) 

As the function |fca(z, ^)p is absolutely integrable, we see immediately that for / 
as in (iii) there is 

limfai(/.)- [ \kaiz,0\''d^zd^^ = ]TT\h^/^-^h'^^] =£a. (6.20) 

We now turn to Ta2 and take into account our assumption 1)6.11) . Using the 
identity (|6.11(l and the fact that 6a and 6ah^^^ are HS, one finds that an equiv- 
alent formulation of the assumption is that h^^^°''> /'^Sa{id +6a)~'^ h^^^"^ ^'^ is a HS 
operator; we denote its kernel in the momentum space by la{—PT(t)- We evaluate 

Ta2{wi,'W2) in momentum space, making use of the identity wT{'p) = uii{—p): 

Ta2{wi,'W2) = j la{p,q)t{p,ct)Wi{p)w2{q)d^pd^q, (6.21) 

where 

m 9) = [{f + + rn-)] (l + -yr^^^^^=^) ■ (6.22) 

As la is square integrable, and t polynomially bounded, 7112 defines a distribution 
Ta2{x,v). Let /,5e 5(R3). Then 

Ta2(e + rj, ^~ ri)f{09irf) d'^ d^ 

= 2 / la{r + s,f-s)t{f+s,r-s)f{2r)g{2s)d^rd^s, (6.23) 



26 



where on the l.h. side the integral notation is symboUc, but on the r.h. side this is 
the ordinary integration. Now, one shows the foUowing estimate 

^(^^+ ^ + 1.1 l^)(3.»)/4 ' (6.24) 

(note that t{p,q) > 0). To prove this it is convenient to consider the cases 
|r|^ > |s|^ + m? and |rp < |s|^ + separately. In the first region one then 
uses the obvious bound t{p, q) < 2[{]y^ + m^){q^ +to^)] "^/^^ while in the second 
one finds that for the given \f\ and |s| the function on the l.h. side is the biggest 
for r-s = 0. Using the bound one easily shows that 

J [t{r + s,f~s)fd^s< const. If]"-". (6.25) 

Therefore t(f+ s, f- s)/(2f) e L'^(R^), so 

la{f+ s, f- s) t{f+ s,f- s)/(2r) e L^{R^) . (6.26) 
Using this fact in 1)6. 23() one finds that 

Sa2{r], f)= I Ta2{^+ V, ^ - 7^) f d^^ 

= -^=i^ J la(r + s,r~ s)t{r + s,r- s)f{2f)e~'^'^^' ^d^rd^s (6.27) 
indeed defines a distribution and is continuous in ff. Thus 

£a2{f) = -ij-^ I lair + s,r-s)t{r + s,r-s}f{2r)d\dh. (6.28) 



\/27r3/2 . 

Using square-integrability of la and the estimate (|6.25(l we have 

|fa2(/)|' < const. /|/(r)|Vl'~"rf'^- (6.29) 



For / as in (iii) one easily then finds 

lim£,2(/e) = 0, (6.30) 

e—>Q 

which ends the proof of our claims. 

In our calculation of the energy density we have used the standard definition 
of the Wick normal ordering. As an aside, it may be of interest to mention that 
this definition has been recently improved for the cases where the reference state f2 
depends on external fields (as e.g. in a fixed curved classical spacetime). The 
problem with the usual definition in such cases is, that the scalar subtraction 
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function depends nonlocally on the background. This may be remedied, as it turns 
out, by an additional subtraction of a smooth function (in "Hadamard states"; see 
the papers by Hollands and Wald ^S], and Brunetti, Fredenhagen and Verch 
for an application to external field electrodynamics see also j^). This has no 
immediate bearing on the discussion in the present work, but may have applications 
in related problems with external fields present from the start (as in a curved 
spacetime) . 

7 Remarks on relations with some other 
approaches 

In this section we make some remarks on the relation of our approach to other 
calculations of Casimir energy. We shall discuss a few characteristic examples of 
local calculations, and next comment on the "zero point" ideology. 

First, to make our point on local Casimir energy, we need to consider a gen- 
eral situation briefly sketched in the Introduction, where two representations of 
local algebras in some open region A^o in spacetime are locally quasiequivalent. 
Suppose we have two representations tt and tt of the algebras of observables in A^o, 
acting in Hilbert spaces TC and H respectively. We assume that the representations 
are locally quasiequivalent, but say nothing on their (global) equivalence. This is 
the expected state of affairs in many situations typically considered for Casimir 
problems. For instance, for a scalar field Mo may be the whole spacetime outside 
some 2-surfaces in 3-space, tt the vacuum representation of the field, and tt the 
representation built on the ground state of the field in presence of the boundary 
conditions imposed on the boundaries of A^o- We choose a state in the representa- 
tion TT, that is a density operator p in 7Y. If the representations are not equivalent 
it makes no sense to ask for a state in the representation tt which gives the same 
expectation values as p for all observables. However, the local quasiequivalence 
tells us that if we restrict attention to an open subset O with a compact closure 
contained in A^o, then there exists a density operator po in Ti. such that 

TT[pn{A)] = TT[pon{A)] for A in O . (7.1) 

The local energy density is not one of the fundamental local observables A, but it 
may be locally built with the use of them. In the sequel we restrict attention to 
the scalar field and construct local energy density as in H6.2(l . Thus given a state p 
the Casimir energy density in O according to the views we follow in this paper is 

Six) = Tv[poH{x)] (7.2) 

(expectation value of a fixed, free field energy density operator). Let p be, for sim- 
plicity, the projection operator onto the unit vector Q. Then using (|6.2|l and H7.1|l 
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we can write for the Casimir energy at a given time t = 0: 

£{x) ^- lim 1(51, [P{x)P{x') + VX{x) ■ VX{x') + Tn^X{x)X{x')]n) 
- {n, [P{x)P{x') + WX{x) ■ WX{x') + m^X{x)X{x')]n)^ , 

where X{u) = $(0, u), P{v) = 0), X{u) = $(0, u), P{v) = 0), and $ and 
$ are operators representing the field under tt and tt respectively. Recall that in 
the case discussed in Section there is H — H, — ft, ^ — $q, and one recovers 
the formula obtained at the beginning of that section. More generally, let in each 
of the representations a different time evolution be given by unitary operators: 
free evolution U{t) and evolution influenced by background U(t) respectively, and 
denote 



Xt{x) = U{t)X{x)U{t)* = ip{t, x) , Pt{x) ^ U{t)P{x)U{t)* 
X^{S) = U{t)X{x)U{ty = ip{t, x) , Pt{x) = U{t)P{x)U{t)* 



(7.4) 



If for both evolutions there is Pt{x) — dLp{t,x)/dt, Pt{x) — dip{t,x)/dt, then one 
can write the last formula for £{x) at i = as 

£{x) = I lim {dtdt'+^-^' + mH{{n,ip{t,x)'f{t',x')n)~{n,ip{t,x)ip{t',x')n)] 
2 t.t'^o 

x' — 

(7.5) 

This formula was derived along similar lines by Kay ITJ in the context of the free 
field in a locally flat spacetime with nontrivial topology. There are no boundaries in 
that case, but the net of local algebras of observables in this spacetime differs from 
that in the globally flat Minkowski spacetime, so the notion of a global Casimir 
energy in the sense we use here has no application. 

In the context of electromagnetic field bounded by conductors in Minkowski 
space the opinion similar to ours, that one should compare expectation values of 
the fixed free field energy density, was expressed by Scharf and Wreszinski 
Consider a massless scalar field analogy of the setting. Then A^g is the spacetime 
region outside boundaries. Put m = in the last formula, use the translational 
symmetry of the two-point functions (which enables the replacement V —^ — V), 
and the wave equation, which both correlation functions satisfy outside the bound- 
aries. This leaves us with 

£(f) = i lim {5t9t--l-942}{(f2,^(i,f)(^(t',f')f^)-(f^,<^(i,^)<^(i',^')f^)}, (7-6) 

X — ^x 

which is the formula used in ^^l- No global energy density may be obtained in 
this way (if not by an ad hoc regularization of the infinities in the density) due to 
the algebraic problems explained earlier. 
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Next, we want to comment on the "Green function" method. In papers follow- 
ing this method one usually states that the (local) Casimir energy is the difference 
between the energy "in the vacuum state with the barriers present and with them 
absent" (see e.g. ^01); with no further explanation on what energy is meant. Stay- 
ing with the scalar field as our example, one then uses with not much comment 
a formula similar to H7.5|l , in which, however, the products of fields are replaced 
by time-ordered products. This brings no change of the result in this simple case, 
but in general has to be justified. As long as outside the barriers the field follows 
the same local equation (the distinct time evolutions agree locally), the ambigu- 
ity as to what energy is meant does not show up. However, this does not matter 
only because for sharp boundaries one cannot determine the global energy any- 
way. And in fact, if the barriers are replaced by external fields one has to make it 
clear what is being calculated. An example of such calculation is attempted in [5], 
where one of the sections treats on the quantum Dirac field in an external classical 
electromagnetic field. The authors' intention apparently is to compare the energy 
of the Dirac field itself, so they keep the free field energy expression. However, they 
take over the form of this expression containing time derivatives, and to eliminate 
them they use different field equations in the two cases, which spoils the original 
intention (remember that the Dirac equation is first order, so the time derivative 
of the field is not an independent initial value variable) . 

Another example of the external field calculations is to be found in . Here 
the authors with the intention of finding the global Casimir energy explicitly com- 
pare expectation values of two different energy operators: energy of the field with 
the interaction terms included and the energy as given by the free field theory, in 
the ground states of the two respective evolutions. In our opinion this is one of the 
reasons for the appearance of infinities in their expressions, which are eliminated 
by adding "counterterms" to the model which does not need them (except for 
trivial normal ordering of quadratic observables). We note, moreover, that it does 
not follow from the smoothness of the external field alone that the ground state 
representations, with the external field present or not, are equivalent globally. 

In the rest of this section we try to understand, from the point of view of the 
formalism presented in this paper, how "zero point" expressions may arise in the 
context of Casimir effect for quantum fields. In our opinion their appearance is 
a consequence of unjustified manipulations. Accordingly, the equations and trans- 
formations to be found below are not to be taken at face value. We indicate this 
by putting a dot over the equality sign. 

The "zero point" expression for Casimir energy has the form 

£l-^-^\Y.i^ak-\Y.u^k. (7.7) 
fc fc 

where uju and ojak are appropriately discretized frequencies of free and perturbed 
field respectively. In our language this would be 

£1-^- ^\TT{K-h), (7.8) 
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which usually is meaningless, but is then "regularized" to squeeze a finite result. 
We show how this expression may arise. 

We have shown in the previous section that £a is a limit of the energy density 
distribution value for the test function tending to one, see Eq. (|6.14|l . Also, it turned 
out that in this hmit only the part £ai{f), Eq. (|6.19|l . of the density distribution 
contributes. Thus we can use Eq. (|6.2l)|) to calculate the total energy. Distribution 
£ai{f) is determined by part Tai, Eq. H6.9|l . of (ila, ■.H2{wi,W2)- ^a)- If we do not 
pay due attention to domains we can rewrite Eq. (|6.9|) by expressing it in terms of 
h and ha instead of h and Sa ■ The result is 

Tai{wi,W2) = j J{x\ha-h + h{h-^-h-^)h\y)wi{x)w2{y)d^xd^y, (7.9) 

which implies 

Sa^^J ^^l'*'^ ^ + ^(^""^ ^ h-^)h\x) dPx 

= '^TT[ha-h + h{h-^ -h-^)h]. (7.10) 

Let us now, again ignoring difhculties, apply this to the case of sharp boundaries, 
where ~ —Ag^^g^-^ with appropriate boundary conditions B{a). Suppose that the 
support of Wi and W2 stays outside the boundaries. Then h'^Wi — h'^Wi, and we 
have 

iW, {ha - h)w2) = \{WI, {h-^ - h-^)h^W2) + ^{m, h^h-^ - h-^)w2) , (7.11) 

(written in two terms only for symmetry reasons), or, for x and y outside the 
boundary, 

{x\ha ~ h\y) = \{x\{ha' - h-')h' + h\h-^ - h-^)\y) . (7.12) 
This needs regularization on the boundaries. Assuming some form of it one writes 

Tr{ha -h) = ^ J{^\{h-' - h-^)h^ + h^h-^ - h-')\x) d^x 

= i Tr [{h-' h-')h' + h^h-' h~')] . (7.13) 

Suppose that the regularization used cuts high momenta, so as to allow one to 
change the order of operators under the trace sign. Then 

Tr{ha-h)^Tr[h{h-^^-h-^)h]. (7.14) 

Using this in H7.10|l one arrives at H7.8|l . 
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A Appendix. Fock space operators and Bogoliu- 
bov transformations 



In the appendix we give a brief review of some known properties of Fock space 
operators which are needed in the main text. The main sources of reference for 
Section lA.ll are books 9: (voLII) and (3). The content of Sections lA. 21 and I A. 31 is 
a rather common knowledge. Precise original proofs of the criterions of equivalence 
of representation use rather more advanced and less common techniques [IE] , so 
we think a simple proof with the use of creation/annihilation operators is worth 
presenting in I A. 41 (The results on equivalence have been later generalized, in the 
widest form in jl9|l. 



A.l Weyl system in a Fock space 

Let TC be the symmetric Fock space based on the "one-particle (excitation) space" 
/C, i.e. 

OO 

W = 0Wn, no = C, = 5(/C . . . /C) (71 > 1) , (A.l) 

""'^ n times 

where S is the symmetrization projection operator. The scalar product in Ti. will 
be denoted by (. , .), the "Fock vacuum" vector by fi, and the particle (excitation) 
number operator by N. On the domain T^i^N^/"^) the annihilation and creation 
operators are defined in the usual way: for each f & K. and -07 X G I?(iV^/^) one 
sets 

a*(/)V = 5(/®V]vTTV), (x,a(/)0) = (a*(/)x,V^), (A.2) 
and shows that 

\\a*{m<\\f\\\\{N+l)'/^^, a#(/)=a(/)ora*(/), (A.3) 
and for ip e V{N) 

[a{f),a*ig)]^^{f,g)ip. (A.4) 

Operators a(/) and a*(/) are respectively antilinear and linear in /. 
Let Tif he the finite-excitation subspace (dense in Ti.), i.e. 

oo k 

U 0^»- 

fc=0 Tl=0 

Operators $o(/) are defined in the following way. One initially sets 

$o(/)V'= ^(a(/) + a*(/))V for ^eV{N'/^). (A.6) 
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Using the bounds ljA.311 one shows that these operators are essentially selfadjoint 

on Hf, so their closures $o(/) are selfadjoint. For -0 € 'D{N'^/'^), ip € 'D{N), 
fi 9, fk ^ K, and real a, /3 one has 

$o(«/ + /3<?)V = a*o(/)V' + /3$o(.9)V' , (A.7) 

a(/)V'-^(<i>o(/) + *<I'o(*/))^, a*(/)^= -^($o(/)-z<i>o(«/))^, (A.8) 

if ^0 then |l$o(/fc)^-$o(/)V'll ^0 (fc ^ oo) , (A.9) 

[Mf),M9)]v^^Mf,9)v- (A.IO) 
Using these relations one shows that the Weyl operators defined by 

Woif) = e^'^oif) (A.ll) 

have the following properties 

W^o(/)M^o(5)-e-5M/.9)H/Q(/ + g)^ VFo(/)* = Wo(-/) , ^^0(0)= id; (A.12) 

the set {Woif) \ f is irreducible; (A.13) 

(r!,T4^o(/)f^) =e~ill^ll'; (A.14) 

if ll/fe-/II^O then ||Wo(/fe)^- Wo(/)^|| ^0 (fc ^ 00) , e H . (A.15) 

Let [/ be a unitary operator in /C. One defines a unitary operator T{U) in Ti. 

by 

r(f/)Wo(/)i^ = Wo{uf)n , (A.f6) 

which implies 

T{U)Wo{f)T{ur - Wo((7/) . (A.f7) 
It is then easy to show that 

T{U) -.Hn^Hn, r([/)|^ ^ U(». . .<»U . (A.18) 

n times 

Let now h he a selfadjoint operator in JC. Then r(e**'*) is a one-parameter group 
of unitary operators. The generator of this group, denoted dr{h), is a selfadjoint 
operator, 

r(e**'') = exp(itdr(/i)) . (A.19) 
Let Vh be any domain of essential selfadjointness of h and denote 

00 k 

Vdvih) = U ^Sj VhQ . . .QVh) . (A.20) 

n tunes 

which means that I'(ir(h) is formed by finite linear combinations of symmetrized 
products of vectors from Vh. One shows that 

dV{h)\gi^^^^ Q^^^ = id® . . . (g) id+ . . . + id® . . . (g) id«)/i, (A. 21) 
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and that dT{h) is essentially selfadjoint on 2?dr(/i)- We assume now that /i is a non- 
negative operator. Then dT{h) is also nonnegative and has the following repre- 
sentation in terms of quadratic forms. As (a(/))* is densely defined (its domain 
contains I?(A^^/^)), the annihilation operator a{f) is closable, we denote its closure 
by d{f). Let {fi} be any orthonormal basis of /C formed of vectors in 'D{h^/'^), and 
denote 

Q{q) = {^A e 7i I V e [\V{d{h^'^U)) and ^ \\-a{h^'^ < 00} . (A.22) 

i i 

One shows that the following form on Q{q) is closed 

X) = Y.^d{h''^h)i^, -a{h''^fr)x) ■ (A.23) 

i 

It is easy to check by direct calculation that the restriction of this form to "Dj^rih) 
gives 

qi^,x) = {^,dTih)x). (A.24) 

As I'(jr(/i) is a core of dT{h), the unique selfadjoint operator defined by the form q 
is identical with dr{h). Thus 

Qiq)^VidT{h)^^^) and q{iP,x) ^ idr{h)^^^ij,dr{h)^/\) . (A.25) 

In particular, for all tp,x^ 'D{dr{h)) identity ljA.24|) holds. 

We note that the particle number operator may be represented as a special 
case of this construction, 

N = dr(id) . (A.26) 

A. 2 Symplectic transformations of (/C, Im(. , .)) 

Hilbert space /C, as a real vector space, is a symplectic space with the form Im(. , .). 
Its real-liner, bijective transformation L is a symplectic transformation if for all 

Im(L/,Lg) = Im(/,g). (A.27) 

The inverse transformation is then also a symplectic transformation satisfying the 
same condition. Substituting / — > L~^f in (|A.27p one has 

Im(/,Lg) = Im(L-V,5). {A.28) 

One defines operators on /C: 

T=\{L~iLi), S^^{L + iLi), L = T + S, (A.29) 
r = ^{L-^ - iL-h) , S' = \{L-^ +iL-^i), L-^=T' + S'. (A.30) 

Operators T and T' are complex-linear, while S and 5" are complex-antilincar. 
Using their definitions and the relation ljA.28|l it is easy to show that operators in 
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the two pairs T', T and S' , -S are mutually adjoint, so T = T** and S = S**. Thus 
both operators are everywhere defined and closed, so they are bounded. Separating 
the identities L^^L — id and LL~^ = id into linear and antilinear parts one gets 

T*T = S*S + id, T*S^S*T, (A.31) 
TT*^SS*+ id, TS* = ST* . (A.32) 

Conversely, if the operators T and S satisfy the above relations on the whole 
Hilbert space JC, they are bounded and define a symplectic transformation 
L = T + S. Furthermore, if the relations are satisfied, then T is a bijection of 
/C onto /C. Thus if T = C/t|T| is its unique polar decomposition, then Ut is a uni- 
tary operator. We set S = UtR- It follows then from the first equalities in HA.31|I 
and (|A.32p that R*R ~ RR*. If R — K\S\ is the unique polar decomposition 
of R, then this condition is equivalent to K\S\ — \S\K, so is a partial anti- 
isometry of (Ker \S\)-^ onto itself. From the first relation in (|A.31|) |Tp = id + |5p, 
so K\T\ = \T\K as well. The second relation in ifOTl) then gives K* = K. 
We summarize the results: 

T ^UT(id+\S\Y^^ , S^Ut\S\K, [\S\,K]=0, (A.33) 
\S\ is bounded, Ut is unitary, if is a conjugation on (KerlS"!)^ . (A. 34) 

Conversely, if these conditions are satisfied, then T and S satisfy conditions ljA.31|l 
and (|A.32p . and determine a symplectic transformation by 

L = T + S, L-^^T*-S*. (A.35) 

If 15*1 has no continuous spectrum, then it follows from the above relations 
that its orthonormal basis of eigenvectors may be chosen such that 

1^1/. = X^f^ , Kf, = /, . (A.36) 



A. 3 Bogoliubov transformations in a Fock space 

With the notation of the foregoing subsections let L = T + 5* be a symplectic 
transformation of the space (/C,Im(., .)), and let us denote 

WoL{f)^Wo{Lf). (A.37) 

It is easy to show that these new operators also satisfy the Weyl relations (|A.12|I . 
The transformation Wo{f) i-^ VFol(/) is called a Bogoliubov transformation. Its 
equivalent form is 

MD^'^odn^MLf), WoL{f)^e'^OLif) _ (A.38) 
For e V{^o{Lf )) n V{^o{Lif )) one defines 

ULim = ^(*ol(/) +^$OL(^/))V', aUf)^ = l={^,^{f) - ^<^>,^{^f))^l; , 

(A.39) 
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and shows by a simple calculation that for ijj £ T){N^/'^): 

aL{f)i^ = a{Tf)^ + a*{Sf)i^, alU)^ = a* {T f)4> + a{S f)^ . (A.40) 
Then using Eqs. lfOT|l one also finds for e V{N'^/'^) 

a{f)^ = aL{T*f)ij-al{S*f)ij, a* {f)ij ^ aUT* f)^: - aL{S* f)i> , (A.41) 
and for Lp e V{N) 

[aL{f),al{g)]ip = {f,g)ip. (A.42) 

One says that the Bogoliubov transformation is implementable in Ti if there 
exists a unitary operator Ul such that either of the following (and then both) 
conditions hold 

WoL{J)^ULWo{f)Ul, <i>oL{f)^UMf)Ul, /e^- (A.43) 

The necessary and sufficient condition for the implementability of the Bo- 
goliubov transformation is that 5 be a Hilbert-Schmidt operator, i.e. 

Tr [S*S] < oo . (A.44) 

If the condition is satisfied, then there exists a unique, up to a phase factor, 
normalized vector Qj^ satisfying the conditions 

aL(/)f^L = 0, felC. (A.45) 

Moreover, one has 

oo 

e f]V{N^/^) . (A.46) 

1=1 

Equations 

ULa*{h)...a*{fk)n = alih)...al{fk)nL, A; = 0,1,..., (A.47) 

with arbitrary test vectors fi, define the unique (up to a phase factor) unitary 
operator Ul implementing the Bogoliubov transformation. For the completeness 
we sketch a simple proof of these statements in the next subsection. 

A slight generalization of the above results is needed in the main text. Let 
J' and J'' be real subspaces of /C, dense in JC, and let L : J' J'' he a. bijective 
symplectic transformation (i.e. a real-linear transformation satisfying Eq. ljA.27p 
for / e j/). Suppose that there exists a unitary operator Ul such that 

WoiLf) = ULWoimi , fej. (A.48) 

Then L and extend to bounded symplectic transformations on JC, and 

Eq. (|A.43p is satisfied for all / € /C. 
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Indeed, suppose that IIA.48P is fulfiUcd. Then using Eq. IIA.14p one finds 

e'lW^fW' = {Uin,Wo{f)Uin), fej, (A.49) 

which shows that L is a continuous transformation on its domain (if /„ ^0, then 
by Eq. (|A.15|) also L/„ 0). Thus L extends by continuity to a bounded operator 
on K. From Eq. (TOHll we have UIWo{I)Ul = Wo(L-V) for / G J', thus similar 
reasoning shows that the extension of L is a bijective symplectic transformation 
of K. onto itself. Equation ljA.48|) now extends by (|A.15|I to all / G /C. 



A. 4 Proof of the statements flA.44l - fA.47p 

Let the Bogoliubov transformation be implemented as in Eq. (jA.43p . For each 
pair of vectors tp,ip <E V{N^^^) one has then {a*^{f)ip,ULv) = {'ip,ULa{f)(p)- 
We substitute here ip = Q, f ~ T^^g, and use Eq. HA.40p . This yields 

{a*ig)^,nL) = ^{a{ST-'g)^,,nL), (A.50) 

where fl^ — Ul^. Substituting here for ip all vectors of the form a*{gi) . . . a* {gk}^ 
for k = 0, 1, . . . recursively, it is easy to see that (f2,f2i) cannot vanish, as oth- 
erwise Hl would be orthogonal to the whole Hilbert space. Let now {fi} be an 
orthonormal basis and put in ljA.50p g — fi and tp = a*{fj)^l, which gives 

{a*{n)a*{f,)n,nL)^--{ST-'f,jj){n,nL). (a.si) 

Take the sum over i,j of the absolute values squared of both sides of this equation. 
On the l.h. side one then gets a quantity smaller or equal 2||r2L|p, so ST~^ is 
a Hilbert-Schmidt operator. As T is a bounded operator, the same is true for S. 

Conversely, let now 5 be a HS operator, so there exists an orthonormal basis 
{fi} satisfying ljA.36|) . and denote gi = Urfi, which defines another orthonormal 
basis. Then 

Sfi = Xtgt, Tf, = Jxf + Igi , so ST^^g,= Z^' g,, 

y' A^ + 1 

(A.52) 

where ^ Af < oo . 

i=l 

We look for a vector Jl^ which lies in the domain of all operators aL{f) and satisfies 
Eq. (|A.45|I . If such vector exists, then it must satisfy Eq. IjA.SOj) for all possible g 
and -0. It is sufficient to substitute for g all basis vectors gi and for ip all vectors 
from the basis of the particle number representation {\ni,n2, . . .)} with profiles 
gi,g2, ■ ■ ■■ This gives the recurrent conditions 



(ni...n, + l...|r!i) = ^L=J-^(ni...n,-l...|r!i)., (A.53) 

V Ai + 
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which are solved for numbers (nin2 . . . \^l) uniquely up to a common constant 
factor c by 

{nin2 . . . ir^i) = if not all are even , 



(2toi 2m2 . . . IQl) = c n f - Jl— ) 



A,; /(2m, -1)!! (A.54) 



Using these explicit expressions one finds that for each I = 0, 1, . . . the following 
sum converges: 



5^ (J2r^^) \{mn2...\nLW 



mi, 7712,... i=l i — l 

The first inequality is obvious, while the second bound will be shown below. For 
I = the bound shows that the coefhcients (nin2 . . . jfii) indeed define a vector fi^ 
solving the conditions ljA.53p . The bounds for I £ N show that this vector is in the 
domain of all operators A^'/^. In particular, fl^ is in the domain of all operators 
o-lU)- This completes the proof of the existence and uniqueness (up to a phase) 
of a normalized vector solving equation (|A.45|I . and of the property l|A.46p . State- 
ments about the operator Ul are now easily proved with the use of commutation 
relations, and the irreducibility of the Weyl system. 

To prove the missing step in Eq. ljA.55|l we denote qi = A?(l + A^)^^. From the 

oo 

finiteness of the sum A^ it follows that also the following expressions converge: 

1=1 



oo 



^'-E(t^) =E^" foraU/GN 

OQ ^ oo / OO 



r — 

i=l ^ 1=1 ^ i=l 

One shows by induction with respect to I that 

^ oc- ^ ; oc 

mi, 7712,... ^7=1 ^ j = l 

where Wi are polynomials. For / = the l.h. side is an infinite product of geometri- 
cal series, so the equality holds with Wq = 1. The step from Z to 1 is obtained by 

°° d 

the application of the homogeneity operator > qi — — to the both sides of the equa- 
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OO ^ OO Q 

tion. A direct calculation yields y qi— — r — pir and > qi— — pk = k(pk +Pk+i), 

^ oqi ^ dqi 

v. ■ i ■ ■ confirms the inductive claim and completes the proof of the bound ljA.55|) . 
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